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We investigate the mechanical stability of a segment of DNA wrapped around a histone in the 
, nucleosome configuration. The assumption underlying this investigation is that the proper model for 

' this packaging arrangement is that of an elastic rod that is free to twist and that writhes subject to 

^SJ I mechanical constraints. We find that the number of constraints required to stabilize the nuclesome 

> I ■ configuration is determined by the length of the segment, the number of times the DNA wraps 

' around the histone spool, and the specific constraints utilized. While it can be shown that four 

constraints suffice, in principle, to insure stability of the nucleosome, a proper choice must be made 
to guarantee the effectiveness of this minimal number. The optimal choice of constraints appears to 
' bear a relation to the existence of a spiral ridge on the surface of the histone octamer. The particular 

04 , configuration that we investigate is related to the 30 nanometer spiral, a higher-order organization 

of DNA in chromatin. 



«4-H 
O 



o 



I 



87.15.By, 87.15.La, 87.15.-v, 62.20.Dc 



I. INTRODUCTION 



' O The issue of DNA packaging has been the subject of intense research for the past forty years. The remarkable fact 

^ that a meter of DNA (the total length of the human genome) fits into a cell nucleus having a typical radius of a 
few microns, and that, thus packed, still manages to perform all its biological functions, has captured the attention 
and continues to challenge the ingenuity of researchers. The fundamental unit of DNA packing in eukaryotes is the 
nucleosome In the nucleosome configuration, a portion of the DNA strand wraps approximately one and three 

■ quarter times around a protein spool, known as a histone [||-^. A string of nucleosomes is believed to participate 
in the next higher order of DNA packing by folding to form the so-called 30- nm fiber Even higher orders of 
' organization have been conjectured, but as yet there is nothing approaching a complete understanding of the physical 
I structure, at all orders, of the DNA in the cell nucleus. Indeed, the detailed oragnization of DNA and proteins in the 
ff^ , 30 nm fiber is not entirely settled ]7|-p^. 

In this paper, our ultimate focus will be on a single nucleosome. We will present a model for the action of histones 
^""^ . in nucleosomes. We treat the segment of DNA in a nucleosome as an elastic rod and apply an approach first developed 
' by Kirchhoff ||ll|] to obtain the equilibrium configurations of DNA in the absence of histones. We assess the stability of 
the elastic rod with respect to small deviations from the equilibrium configurations, and we find that all configurations 
that are not equivalent to a straight rod are unstable to fluctuations. More specifically, all non-trivial equilibrium 
configurations represent saddle-points in energy space. Then we prove that the presence of histones provides a physical 
mechanism by which more compact configurations of DNA are rendered stable against purely mechanical fluctuations. 
(~| The specific mechanism is a set of constraints on the fluctuations about mechanical equilibrium, which can be simply 
O 1 modeled mathematically. The stability of the equilibrium configuration is then framed in terms of the determinant 
O ' of an n-hy-n matrix, where n is the number of constraints. Of special interest to us is the question of the most 
^ , economical combination of constraints that serves to stabilize a given configuration. "Economy" in this case refers 
k>( " number of constraints required to accomplish stabilization. We are able to establish in the case of particular 

rN I interest to us that a nucleosome-like configuration of DNA is rendered stable by four constraints, and that no fewer 

. constraints will accomplish this. 
" " ' The conjecture that underlies the work reported here is that the portion of DNA that is wrapped around a histone 
is in a state of unstable mechanical equilibrium that is rendered mechanically stable by the constraints associated with 
the histone spool. There are reasons to believe that such a state is desirable. Imagine a pencil standing vertically on 
its point. While this state will not persist in the absence of outside influences, it can be sustained without application 
of substantial external forces. In fact, no force at all is required to guarantee the persistence of this state if the pencil 
is exactly vertical. By the same token, the removal of the constraint that keeps the vertically balanced pencil from 
toppling is accomplished at no energetic cost. Viewed in this way, the nucleosome configuration represents a highly 
efflcient stratagem for the local packing of DNA, in that the histone spool is introduced and removed with minimal 
expenditure of resources. 



1 



Of course, because DNA is, in its "naked" state, highly charged, electrostatic interactions play an important role in 
the behavior of this molecule, both in vitro and in vivo. These interactions are, apparently, key to the condensation 
of DNA in prokaryotes Furthermore, electrostatic interactions, highly screened though they may be, could be 
important contributors to the energetics of the DNA-histone interactions. They are ignored in this work. Nevertheless, 
the elucidation of the stabilization of configurations within the context of the purely clastic model of DNA ought to 
aid in the investigation of similar questions for DNA subject to the influence of additional interactions and degrees of 
freedom. 

An outline of the paper is as follows. In Section || the Hamiltonian governing configurations of DNA, modeled as a 
bent and twisted rod, is presented, along with the quadrature solutions to the extremum equation for that Hamiltonian. 
Section [II reviews the formulation of linear stability analysis in this case. A key result of the considerations outlined 
in this section is that any non-trivial extremal configuration of bent and twisted DNA will be mechanically unstable 
if the strand is long enough, and if there are no constraints on fluctuations. This means that all extremum solutions 
for an infinite strand of bent and twisted DNA represent saddle points of the energy. Section IV introduces the 
notion of mechanical constraints, particularly those constraints associated with the requirement that the straight-line 
distance between two points on the DNA does not change as the segment fluctuates. The way in which this and other 
constraints are mathematically implemented is discussed in this section, and in Appendix ^ In section ^ we address 
the issue of the number of constraints required to stabilize a section of bent and twisted rod against fluctuations about 
an extremum solution. We find that the minimum number of constraints needed to do this is equal to the number of 
unstable eigenmodes of a fluctuation operator introduced in Section HI. Section VI contains a discussion of the effects 
of a periodic array of constraints on the stability of an infinitely long section of bent and twisted rod. This discussion 
provides a lead-in to our investigation into the influences required to stabilize the nucleosomal configuration of DNA. 
It is also relevant to the stabilizing action of a protein armature on DNA in chromatin. Sections VII and VIII directly 
address the issue of the stabilization of DNA in the nucleosomal configuration. There are four unstable modes in the 
case that we investigate. We find that four constraints suffice to counteract them. However, we also find that those 
constraints must be chosen with care. The effective set reflects the known structure of histone, in particular a spiral 
groove that has been identified. In our model, this groove acts to limit the ability of the DNA wrapped around it to 
slide parallel to the spool's axis. 



II. HAMILTONIAN AND EXTREMUM SOLUTIONS 



The configuration of the twisted, writhing rod is characterized in terms of the Euler angles, depicted in Figure 
We assume an isotropic rod, characterized by a bending modulus A and a torsional modulus C. The elastic energy 
of the rod, in terms of the Euler angles 9{s), 0(s) and ipis), is given by 
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Here, s is the arclength along the rod. A way to determine the equilibrium configurations of the rod is to supplement 
-^clastic with the term 



Ec^-F J cose{s)ds 



(2.2) 



This contribution can either be seen as a Lagrange multiplier that enforces a given end-to-end distance, or as repre- 
senting the effect of tension on the rod. 

Finally, in certain cases an additional constraint guarantees constancy of the end-to-end linking number. In the 
case of a rod with clamped ends, this quantity is given by 
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Note that (2^) is an integral over a perfect differential. This reflects the topological character of the linking number. 
A fixed linking number is enforced with the use of a Lagrange multiplier. The quantity to be minimized is, then, the 
combination 
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(2.4) 



As the linking number is not a quantity of interest here, we ignore the final term in (2.4). 
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The extremum equations have been extensively investigated They are identical to the equations for the 

behavior of the heavy symmetric top. The connection between top motion and the deformations of a thin rod was 



first noted by Kirchhoff, in 1859 |11|, and is known as the "kinetic analogue" Those equations reduce to the 

following set of three: 

d(j) _ - J^, cos 9 ^2 g-j 



ds A sin^ 9 

cose* (2.6) 




ds C ds 

Fcos9 (2.7) 

The quantities J^, and Eq in the above equations are integration constants. The quadrature result for the angle 
9 is the result of integration of SEtot/S9. Defining u = cos 9, the behavior of solutions can be extracted from the 
equation for 9: 

ds = ^ (2.8) 

(^Eo -Fu)-i, ( J2 + Jl^ 2J^J^u) y/'-^{u-a){u-b){u- c) 

where c < u < b < a. The characteristics of the solutions depend on quantities a, b and c. The property of solutions 
have been completely investigated ||l5|| for different values of a, b and c. We focus our attention on those configurations 
which have the same configurational form as a segment of DNA in a nucleosome. For specific value of a, b and c we 
have obtained a solution which is depicted in Figure ^j. The reason that we have made the choice of parameters above 
has to do with the close visual relationship between the conformation of the bend and twisted rod as displayed in Fig. 
^ and the DNA in a commonly conjectured form of the 30-nm spiral 1^. We operate here under the assumption that 
this configuration is a reasonable representation of the organization of DNA in this component of chromatin. 



III. STABILITY DETERMINATION 



The stability of an extremum configuration is determined by altering the configuration and calculating the change 
in the quantity that is extremized. In this case, the quantity of interest is the elastic energy. A stable solution is one 
that minimizes the energy. If the quadratic effect of any small deviation from this solution is to lower the energy, 
then the solution cannot represent stable equilibrium. Instead, the configuration is unstable; it is either a maximum 
energy configuration, or a configuration at a saddle point of the energy. In the case of the equations for the elastic 
energy of a twisted and bent rod, the second order effect of second order fluctuations is obtained by taking second 
functional derivatives of the expression in Eq. (|]^) with respect to the Euler angles 9, (j) and ip, and then by setting 



those Euler angles equal to their classical values. After a bit of reduction, we find that the question of the stability 
of a classical configuration can be framed in terms of the spectrum of the following operator 

Vis) (3.1) 



where 

1 x{s) 1 
4(l-u(s)2)2 ^ 2 

and 



x{s) = (2 - u{s)) (1 - u{s)f [a + l){b + l)(c + 1) 

+ (2 + u{s)) (1 + u{s) f [a - 1)(1 - b){l - c) (3.3) 



Here, u{s) is the solution for cos 6'(s) displayed implicitly in Eq. ( |2.8[ ). If all eigenvalues of the operator (3.1) are 
positive, then the solution to the classical equation is stable. If any is negative, then there are fiuctuations that 
decrease the energy to a value below its classical value. Note that the operator in question resembles the Hamiltonian 
for a one-dimensional particle in the potential V{s). In figure ||, we display this potential for the choice of parameters 
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a, b and c utilized in this investigation. The potential, which has been displayed for an extended range of the arclength, 
s, has the form of the sort of periodic potential encountered in discussions of electrons in metals and semiconductors. 
As in this case, the eigenvalue spectrum consists of bands of "allowed" states, separated by "forbidden regions." The 
question is whether all or part of any of the bands lie below zero. As it turns out, this is, indeed, the case. 

The reason for this is that one can identify a mode having an eigenvalue that is strictly equal to zero — the trans- 
lational mode, equal to the derivative with respect to s of the classical solution, dc\{s). The existence of this mode 
follows from the translational invariance of the extremum equations, and it is known to play a key role in, for instance, 
the question of tunneling between the false and the true vacuum in quantum field theories |17|| . The translational 
mode in this case is displayed in Figure |^. Note that this mode is not spatially uniform, and, in particular, that 
it possesses nodes. On the basis of elementary considerations, one knows that there result, solutions to 

the effective Schrodinger's equation associated with lower — hence negative — eigenvalues. Figure || displays the band 
structure associated with the potential in Figure ^. In an infinitely long section of twisted and bent rod that has taken 
the configuration pictured in Figure ^ there is an infinitely large set of distortions that will lead to a lowering of the 
rod's total elastic energy. In fact, the reasonable expectation is that these distortions are a route to interwinding. 



IV. THE SOURCES AND MATHEMATICAL IMPLEMENTATION OF PHYSICAL CONSTRAINTS 

Stabilization of the classical, or extremum, configuration can be achieved by the introduction of mechanical con- 
straints. These constraints are expressed mathematically as the requirement that a property of the DNA's configura- 
tion does not change under distortions about the extremum solution. The physical constraint that a certain quantity 
be kept constant translates fairly straightforwardly into a set of mathematical conditions on the fiuctuation spectrum. 
In turn, these mathematical restrictions lead to a reformulation of the method by which the fluctuation spectrum 
is determined. The reasoning leading from physical constraints to a new approach to the determination of effective 
eigenvalues of the linear fiuctuation energy operator is presented in Appendix ^ 

Briefiy, a constraint on the conformation of a bent and twisted rod is expressed mathematically in terms of a 
condition of the form 

T {0{s), (j){s),4'{s)) ds = a constant (4-1) 



Under the assumption that Eq. (4T) is satisfied for the extremal configuration, one then expands to the first order 
in deviations from the extremal forms of 0{s), 0(s) and ipi-s). As it turns out, the first order corrections to (j){s) and 
ip^s) are readily expressed in terms of the correction to 0{s). This is because of the simple way in which these two 



angles enter into the expression for the total energy i n E q. (2.1). If we denote by 7(5) the displacement of 9{s) from 



its "classical" form, the general constraint equation, (4T) becomes 



f{Oci{s), Ms)h{s)ds = (4.2) 

where the subscript "cl" indicates that the quantity is a solution to the extremum equations. 

For closed configurations, such constraints arise naturally from the preservation of the topology (i.e. linking number) 
of the original configuration, and from the requirement that the distorted rod continues to close smoothly on itself. 
Here, such considerations do not necessarily apply, although one might imagine cases in which a "pinning" of the 
ends of a segment forbids any alteration of the linking number. Nevertheless, "boundary conditions" that result from 
physical constraints on the end-points of a given segment do give rise to mathematical constraints on fiuctuations 
about a given configuration. Whether or not boundary conditions will stabilize a segment of bent and twisted rod 
against thermally-driven fluctuations depends on the length of the segment. If the segment is short enough compared 
to the persistence length of the rod, such stabilization is possible. 

Constraints may also be imposed as the result of physical barriers. For example, imagine that the displacement 
vector between two points on the bent and twisted rod is not allowed to vary. One might imagine such a constraint 
being enforced with the use of a stiff, inextensible "brace" flrmly attached to the rod at the two points in question. 
This brace is then immobilized against rotations. The displacement vector between the two points is represented as 

i?o = xxo + yyo + zzq (4.3) 

where 
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/■S2 

xo = sin 6 {s) COS (j){s)ds (4.4) 

Jsi 

/S2 

sm9{s) sin (f){s)ds (4.5) 
_ 1 

zo= COS 6»(s)rfs (4.6) 

J Si 

The constancy of each component of this displacement vector is ensured by a set of three constraints on the deviations 
of the Euler angle 9{s) from its extremum value. If we write 

eis) = Oeiis) = j{s) (4.7) 

then, the following three conditions hold: 

1. Constancy of the a;-component of the displacement vector 

S.= r Ls)cosMs) + J^ '' ' "!f + "^•^•^^' ^.(^)| ^is)ds = O (4.8) 

J Si \ ZA (1 _ u{sy) ' J 

where the quantity Ix(s) is given by 

Ms)= [ Vl-u{s')^sin(Pciis')ds' (4.9) 

2. Constancy of the ^/-component of the displacement vector 



X I ) ( \ ■ J. t \ rf^Pl(l-"(s))^+P2(l+u(s))^ 1 , \J n (A^C^\ 

Sy= I <; u(s)sm(/)ei(s) - y — (i _ n(s)2)3/2 Iy{s) '> -/{s) ds = (4.10) 



where 



Ty{s)= f ^/l -u{s'f cos (Palis') ds' (4.11) 

J Sl 

3. constancy of the ^;-component of the displacement vector 

V'l - u(s)27(s) ds = (4.12) 

_ 1 

In the above relations 

P(.)^[(c±l)(6±l)(a±l)]^/^ (4.13) 

As an alternative to the above set of three constraints, one might imagine that the projection of the displacement 
vector in a given direction is held constant, in which case the contstraint is a linear combination of those constraints: 

xo6x + yo5y + zqSz = (4-14) 

As an example of the use of this less restrictive constraint on fluctuations of the bent and twisted rod, imagine that 
the brace is allowed to rotate, but that it remains stiff and inextensible. Then the single constraint that holds is 
that the projection of the displacement vector along the original direction of the brace is held fixed. Thus, physical 
constraints on the possible contortions of a strand of DNA translate straightforwardly into mathematical constraints 
on the fluctuations of that strand about its "classical" configurations. 
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V. NUMBER OF CONSTRAINTS REQUIRED TO STABILIZE A CONFIGURATION 



It appears intuitively obvious that the stabilization of a given configuration, when there are a given number of 
unstable modes, requires the imposition of an equal number of constraints. It is fairly straightforwardly demonstrated 
that if there are n unstable modes, then at least n constraints are required to stabilize them. To see that this is true, 
suppose that the operator L has four negative eigenvalues. Also, imagine that three constraints have been imposed, 
of the form 



{f\x,)= J f{x)xA^)dx^Q (5.1) 

for any fluctuation, f{x). Here 1 < J < 3. The four eigenfunctions having negative eigenvalues will be ^i(x), with 
1 < i < 4. Let's define 



9r3 = J (ti^)Xj{x)dx (5.2) 

Now, take a fluctuation that is of the form 

4 

/(a;)=^a,6(a:) (5.3) 

1=1 

The three constraint equations are of the form 

4 

"■^9^j = (5.4) 

i=l 

These are three equations in the four unknowns a^. We'll assume that not all g^'s are equal to zero for any i. Then, 
it is possible to set one of the a^'s equal to one. The equations reduce to three linear, inhomogeneous, equations in 
three unknowns. Unless there is some degeneracy, it will be possible to find a solution to those equations. This means 



that a function of the form (5.3) will obey the constraints. Furthermore the expectation value 

{f\L\f) = J f(x)L{x,x')f{x')dx dx' (5.5) 

will be given by 

4 

(/|L|/)=^a?A. (5.6) 

1=1 

Given that the four A^'s in the sum are all negative, we have a fiuctuation for which the expectation value of the linear 
operator L is negative. 

It is, thus, clear that three constraints do not suffice to stabilize a classical configuration against fiuctuation when 
there are four unstable modes. This conclusion generalizes straightforwardly to the case of n unstable modes and 
m < n constraints. On the other hand, n constraints may or may not prove sufficient to guarantee stability. Consider, 
for example the case of a single instability. Let the unstable mode be £,o{x). If the single constraint requires that all 
fluctuations be orthogonal to x{^)i then the equation satisfied by the eigenvalues, A of the constrained fluctuation 
operator is 



E 



aT^a - ° ^^-^^ 



It is straightforwardly verified that solutions of this equation lie between consecutive eigenvalues, Xi of the uncon- 
strained fluctuation operator. Thus, the lowest allowed eigenvalue of the constrained operator lies above the lowest 
eigenalue in the unconstrained system. However, it also lies below the next-lowest unconstrained-system eigenvalue. 

To see that stabilization may or may not occur in this case, we consider two particular subsets of the many possible 
alternatives for the function xi^)- First, imagine that x(a;) cx ^o(a^)- Then the constraint entirely eliminates the 
unstable mode and stability is guaranteed. On the other hand, suppose that xi^) ^ii^) with i ^ 0. Then, a stable 
mode is eliminated, and the unrestrained unstable mode contributes to the fluctuation spectrum. The instability is 
entirely unaffected. 
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VI. THE EFFECTS OF A PERIODIC ARRAY OF CONSTRAINTS 



In the case of DNA confined to the nucleus of a cell, it is widely conjectured that the packing of DNA is accomplished 
with the use of a hierarchical organization of the long strands that constitutes the genome. Given that this organization 
will not be a mechanically stable structure, at least within the bent-and-twisted rod model of DNA, some constraining 
mechanism is required. The histone spools provide one set of constraints, but these operate on the lowest level of 
organization. It is possible that a protein "armature" provides the necessary stabilization at higher levels. Here, we 
discuss the implications of the kinds of mathematical constraints on fluctuations about unstable mechanical equilibrium 
that one can reasonably associate with the mechanical influence of this mechanism for stabilization. 

In this context, we focus on the case of a long strand of distorted DNA, or, equivalently, a long section of bent and 
twisted rod. Here, the set of fluctuations that lowers the energy of the unstabilized configuration is quite large. This 
implies the need for a large number of constraints. When the strand is infinitely long, and the number of unstable 
mechanical modes is infinite, then an infinite number of constraints is required. We will look here at the stabilizing 
effect of a periodic array of constraints. 

The operator L controlling the stability of the equilibrium configuration of a long segment of the bent and twisted 
rod has the form 

-^+Vi.) (6.1) 

Where the potential term, V{x), is periodic, in that 

V{x + a) = V{x) (6.2) 

According to Floquet's theorem, the eigenfunctions of the above operator are of the form 

Kk{x) = e"="C„,fe(x) (6.3) 

where k, called the crystal momentum in solid state physics, is confined to a Brillouin zone. The most convenient 
Brillouin zone for our purposes is —n/a < k < it /a. The function Cn{x) is periodic in x, in that 

Cn,k{^ + a) = Cn,k{x) (6.4) 

The integer n is called the band index. The eigenvalue of this eigenfunction is also indexed by the crystal momentum 
and the band index, i.e. En,k- 

If a periodic array of constraints is imposed, in that we require all fluctuations to be equal to zero at x = & + ma, 
with m an integer and — oo < m < oo, then the requirement that the determinant is equal to zero translates into the 
requirement that the following product is equal to zero: 

7r/a 

n ^(fc) (6-5) 

fc=— 7r/a 

where 

Hk) = E (6-6) 

^ En,k - A 

This tells us that for every value of fc, the lowest value of A corresponding to a fluctuation lies between the lowest 
value of En^k, as a function of the band index n and the next lowest value of that eigenvalue. Suppose we impose two 
constraints, by requiring that the fluctuations are zero at x = 6i + ma and x = 62 + ma. Then, the determinant will 
consist of a product of terms of the form 



E 



(gn„.A:(?>l)C„,.,(?>2)-C,„.,i?>2)C„^A:(?>l))' ^g_^^ 
na>n, (^"-'^ ~ ^) ^^^^^^ - ^) 

In this case, it is possible that the lowest solution of the characteristic equation will lie even higher than when there 
is only one constraint per period. 

As an indication of the effect of an array of constraints, we consider the case of the eigenstates of the "periodic" 
potential that is equal to zero everywhere. As is well-known, one can imagine a one-dimensional Brillouin zone of fixed 
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width. The dispersion relation can then be expressed in terms of a series of curves in which the "crystal momentum" 
is restricted to this zone. There are no band gaps, but otherwise the bands are well-behaved. Now, one is interested 
in the expectation values of the operator 



(6.8) 



Imagine that the lattice spacing is one, and take for the function to which fluctuations are orthogonal a gaussian of 
the form 



(6.9) 



The array of functions are copies of ( |6.9| ) centered about the points a; = ±1, ±2, ±3 . . . The solution to the equation 
setting (6.5) equal to zero is graphed in Figure ^j. Also shown in that figure are the "bands" of the unconstrained 
operator. Note that the allowed values of A for fixed crystal momentum, g, lie between successive bands. This is a 
general feature of array of constraints. 

The effects of a pair of constraints in every period is illustrated in Figure |^, where the values of the A for both one 
and two constraints per period are compared with the A-versus-q relationship for the unconstrained operator. 

This brings us to the way in which a physical armature, in the form of a protein scaffolding, can act to stabilize 
a nontrivially supercoiled DNA configuration. We imagine a configuration as depicted in Figure ^. The contacts 
between the DNA and the armature will stabilize the DNA against fluctuations. 



VII. THE CASE OF A NUCLEOSOME: PRELIMINARIES 

In the nucleosome configuration a segment of DNA wraps around a collection of proteins known as a histone . 
In the schematic depiction of the nucleosome, the DNA segment is represented as a spiral surrounding a cylinder. See 
Figure ||. As a first step in our investigation of the stability of the spatial configuration of the segment of DNA that 
participates in the nucleosome we will look at the stability of the spiral solution to the energy extremum equations 
for a bent and twisted rod. 

Now, the spiral is a special, limiting case of the solutions to the classical equation for u{s) — cos 9{s). In this 
solution, the equality b = c holds, and u{s), which lies between those two parameters in the classical solution is, thus, 
fixed at their common value, which we henceforth will call b. The complete determination of the solution requires 



that we set the parameter a and choose signs in (4.13). We find that there are four possibilities for the solution, 
corresponding to the four choices of the two signs. Two of the solutions are for a left-handed spiral, and in the 
other two the spiral is right-handed. In the case of the nucleosome, DNA is wrapped around the histone spool in 
a left-handed spiral. Given the sense of the helical solution, there the two alternative solutions are spirals that the 
arclength of a single turn of which is either greater or less than 2Tr y/F/A. The quantity y/F/A is the persistence 
length of the rod, and the only intrinsic length scale in this system. If we rescale arclengths so that they are expressed 
in units of ^JFjA, then the rate of change of the Euler angle 0(s) in this classical solution is given by 



d<i>{s) _ Va + 1 ± Va ^ 



ds 

The diameter of the cylindrical region encircled by the helical solution is given by 

2V2 



- , sin6>= V 1-62^2 (\/^nT\/^^) (7.2) 



while the distance between successive turns of the helix, measured along the direction parallel to the cylinder's axis, 
is given by 

^ 008 61 = V27^(^^TTTV^^)^ (7.3) 

V a + 1 ± V a - 1 

The linear equation, the eigenvalues of which yield the energies of fluctuations about the classical solution, is 

d'^^(s) 

-^^ + {b-a)m^ms) (7.4) 
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The parameter a must be greater than one to ensure real solutions to the classical equations, while h = cob 9 lies 
between 1 and -1. If the length of the spir al is allowed to become infinite, then, in the absence of constraints, there 
are an infinite number of solutions to Eq. (7.4) with negative values of A. 

Eq. (7.4) is just the kind of equation for which one might envision stabilization as the result of the imposition of 
a regular array of constraints. Here, we ask what physical constraints will have the effect of stabilizing the extended 
helix against fluctuations. One possibility is depicted in Figure |o[ The vertical dark lines in the figure represent rods 
that enforce a fixed spacing between a point on the spiral and the point immediately above or below it. As shown in 
the figure, there are two "lines" of these rods, on opposite sides of the spiral. As we will see, this arrangement proves 
sufficient to stabilize a family of spirals against fluctuations. 

The requirement that the distance between a point on the spiral and a point separated from it by a single turn of 
the spiral translates into the following requirement on a fluctuation, <i>(s) 



(l){s)ds = 



(7.5) 



Here, sq is the location of the first point along the spiral, while r is the "period" of the spiral, the backbone distance 
from a point on it to the point one turn of the spiral subsequent. In this case 



V2 



Va + 1 ± Va - 1 



27r 



(7.6) 



We will henceforth take the sign in (7.6) to be the upper one, corresponding to the more tightly wound of the two 
branches. Then, the periodic set of constraints indicated in Figure pleads to the following equation for the eigenvalues 
of the constrained spiral 



E 



1 



where 



{k + 2nujf [k + 2nujf - {X + a - b) 



2n 



= 



(7.7) 



This equation is a specific realization of ( |6.6D , in which factors that are independent of the summation variable n have 
been ommitted. The sum in (7.7) can be performed with the use of contour integration. The equation that results is 



TT W COt(|^ 



W|E5+A) _ cot(|^ + --/™ ) - Va-fe + Acsc(|f) 



= 



4(a-fe + A)^cj2 

The minimum value of the A that solves this equation corresponds to k — iw, at which point 



(7.9) 



A 



Va + 1 + Va - 1 



(7.10) 



A graph of (A + a — h)/uj^ as a function of k/uj is shown in Figure |ll|. The constraints depicted in Figure |l^ will keep 
the spiral in place against thermal fluctuations. 



VIII. STABILIZATION OF A SINGLE NUCLEOSOME IN A 30 NM-SPIRAL-LIKE ARRAY 

Here, we take the point of view that there is merit to the notion of an organized and orderly array of nucleosomes 
in the 30 nanometer spiral, and we search for this order in the solution to the energy extremum equations for a bent 
and twisted rod. Interestingly, solutions that mimic a conjectured form of this higher order structure can be found. 
One such solution is depicted in Figure ||. As previously noted, this solution bears a visual relationship to the coiling 
of DNA in a conjectured form of the higher order structure known as the 30 nm spiral. The speciflc values of the 
parameters a, h and c that generate this configuration are 
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a=l (8.1) 
b = 0.656009130822 (8.2) 
c=-0.85 (8.3) 

In this study, we focus on a particular portion of this structure, corresponding to two loops in the distorted DNA 
strand. This portion is illustrated in Figure |lj. Note that the loops are not compact as in the standard picture of a 
nucleosome. The case here is a bit figurative, as we are interested in the notion of organization on a larger scale as 
envisioned in some versions of the 30 nm spiral. 

The "bare" stability of the two-loop portion of DNA was calculated by assuming that fluctuations were consistent 
with free boundary conditions, in which the slope of the fluctuations in the angle, 6{s) is set equal to zero at the 
two ends of the DNA segment. With these boundary conditions, we find that there are four unstable modes of the 



fluctuation operator (BJJ), with potential V{s) as given by (3.2). The eigenfunctions associated with those fluctuations 



are shown in Figure 14l In line with the discussion above, this implies the need for at least four constraints on 
fluctuations of the segment of DNA that is wrapped about the histone in this configuration. The construction of 
these eigenfunctions required an elaboration of the integration method that we generally utilized to find the solution 
of the linear second order equation that governs fluctuations about extremal solutions. This elaboration is discussed 
in Appendix 

We choose to assume that the histone provides constraints in the most "efficient" manner, that is, that number of 
constraints that follow from the presence of the histone does not exceed the minimum number required to guarantee 
stability of the nucleosome configuration. Histones keep the two loops close to each other and limit the arbitrary 
fluctuations of two loops with respect to each other. With this in mind, we started by fixing the distance between two 
different points on the segment of DNA which wraps around the histone octomer. As shown in the previous section, 
at least four constraint functions are required to stabilize the nucleosome structure. 

Our strategy is to construct four constraint functions, each associated with fixing a different distance on a segment 
of DNA in Nucleosomes. We are then faced with t he p roblem of solving for zeros of the determinant of the matrix 
Gki- This matrix is defined in Appendix in Eq. ( |A6| ). The operator C, is given here by 

^ ' — nT) — ^^■^> 



Here, $(s) is an eigenvalue of the operator (3.1) that has the property <i>'(0) = 0. The quantity L is the total arclength 
of the nucleosomal segment. 

As an initial attempt, we fixed four "diagonal" distances between two loops. We fixed these distances only in x-y 
plane. We assume that histone has a distorted cylindrical shape and this way we fixed the radius of cylinder in four 
different places. With this set of constraints, the DNA segment has some freedom to move vertically as long as it is 
wrapped around the histone spool. In this case, we found that constraints only removed two negative eigenvalues and 
the system remains mechanically unstable. We then tried quite a few set of constraints related to keemng the segment 
of DNA loosely on the histone octomer. A few of these sets of constraints can be seen in the Figure |lj. None of these 
sets of constraints was able to eliminate all negative eigenvalues. An example of the determinant |Gfci|, defined in 
Appendix associated with one of the sets of inadequate constraints is shown in Figure |l6|. 

As indicated by the brief account above, the task of constructing such constraints is by no means trivial. Four 
constraints chosen at random, will not, in our experience prove adequate to the task of stabilizing the DNA segment 
against fiuctuations. To understand the mechanism of removing of a negative eigenvalue better, we constructed a five 
by five matrix with the same eigenvalues as the five lowest eigenvalues of our problem. We let the computer choose 
four constraints randomly and ran the program many times. We were not able to see even one case in which the 
constraints remove the four negative eigenvalues. The distance between the third eigenvalue and fourth (as shown in 
the picture) is very large compared to distance between other eigenvalues. As a result it is not at all easy to find a 
set of constraints that eliminates all negative eigenvalues. 

In the end, consideration of the detailed structure of the nucleosome, and a knowledge of the nature of the periodic 
constraints that stabilize a long spiral of DNA led to four constraints that give rise to mechanical stability |^ . Three 
of the four constraints corresponded to rods that stabilize the segment against motion parallel to the (curved) axis of 
the histone, and the fourth is in the form of a "diagonal" strut, reaching nearly across the double-looped segment. In 
more detail, the constraints correspond to rigid, but hinged, rods that join points in the nuclcosome-like segment as 
follows: 

1. The diagonal strut reaches from a quarter of the way in the first loop, to a quarter of the way from the end 
of the second loop. This is the long, diagonal support illustrated in Figure [l7|. Note that the picture of the 
nucleosome here is figurative, in that the "real" histone, as shown in Figure U3, has a curved axis, so as to fit 
into the loops of the nucleosomal DNA. 
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2. The second support runs parallel to the axis of the histone, from a quarter in the first loop to a quarter in the 
second. This is the topmost horizontal support in Figure [l^. 

3. The third support, also parallel to the axis of the histone extends from the halfway point of the first loop to the 
halfway point of the second loop. 

4. Finally, the fourth support, which, like the second and third ones, runs parallel to the histone's axis, joins the 
point three quarters of the way into the first loop to the point a quarter of the way into the second loop from 
the opposite end. This is the bottom horizontal support in Figure 

The equation for eigenvalues of the fluctuation spectrum, now has the form of the characteristic equation of the 
appropriate version of the matrix Gki displayed in Eq. ( [A6| ). The determinant of this matrix, as a function of 
the eigenvalue parameter. A, is shown in Figure The zeros of the determinant occur at the eigenvalues of the 
constrained fluctuation spectrum. We note that there are no negative roots for negative values of A. The poles that 
appear in the plot lie at the locations of the eigenvalues of the unconstrained spectrum. The four negative energies are 
readily identified in the figure. It is worth noting that the stabilization leaves the segment with a positive eigenvalue 
that lies close to zero. In other words the four constraints that were utilized were adequate to achieve mechanical 
stability, but only barely so. 

As noted above, the choice of the four constraints that led to stability of the nucleosomal configuration was guided 
by known properties of the histone octamer. A variety of investigations has revealed the existence of a spiral "trough" 
in the surface of the histone p| p9|p0[ | . Such a trough will act to constrain wrapped DNA against movement along the 
surface of the histone spool that is parallel to that spool's axis. In particular, the section of DNA that is wrapped 
about the histone spool will not be allowed to move in such a way as to alter the distance between adjacent coils, 
when that distance is measured along a direction parallel to the spool axis. In addition, we were guided by the results 



reported in Section VII, in which it was demonstrated that an extended spiral is stabilized by a periodic array of 
constraints equivalent to a set of rigid, but hinged, rods running separating consecutive turns of the spiral as indicated 
in Figure |l^. 

From investigations of the chemical electrostatic and conformational structure of the histone octamer, it is clear 
that the points of contact between the histone and the DNA wrapped around it exceed the minimal number that, 
according to our results, stabilize the DNA segment against mechanical instabilities. However, it is satisfying that 
a "minimal" set of constraints will also do the job. The significance of this result for the mechanics and biology of 
the nuclesome configuration remains to be worked out. Nevertheless, it has long been known that a few points of 
contact between DNA and the histone spool suffice to stabilize the nucleosome [g2| . We believe that issues of optimal 
efhciency will prove relevant in discussions of the nuclesome in eukaryotic chromatin. 
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APPENDIX A: THE INFLUENCE OF CONSTRAINTS ON THE FLUCTUATION SPECTRUM: 

GENERAL RESULTS 



The mathematical effects of constraints on the fluctuation spectrum of the operator (3.1) are readily expressed in 
terms of the roots of a determinant. Here, we outline the way in which this formulation of the stability investigation 
is arrived at. The discussion in this section has appeared before ||T^. It is repeated here for the convenience of the 
reader. 

The investigation of the stability of a solution to an Euler-Lagrange equation, such as the one relevant to the 
configurations of interest to us here can be framed in terms of the eigenvalue spectrum of a linear operator. This, in 
turn, can be recast in terms of the problem of finding extremal values for the expectation value 



mo 



(Al) 



where C is the linear operator. In the case at hand, L is the operator in ( p.!] ). The constraints are equivalent to 
requiring that the ^ between which the operator is sandwiched is orthogonal to a set of m x's. There is also the 
constraint on the absolute magnitude of ^. The constraints are, then of the form 

m) = 1 (A2) 

(ClxO = (A3) 

In Eq. (A3), the index I runs from 1 to m. The equation for the extremum of the quadratic form ( |Al[ ), subject to 
the constraints (A2) and (A3), takes the form 



1=1 



AilxJ 



(A4) 



The coefhcients A and A; are Lagrange multipliers, which enforce the constraints to which the system is subject. The 
solution to the above equation is 



10 = 



A, 



ti^ 



A 



\Xi) 



(A5) 



The Lagrange multipliers A; must now be adjusted to ensure the orthogonality requirements. These requirements are 
of the form 



= EA^(X.I^ 

1=1 

= GkiAi 



Xi) 



(A6) 



This set of m equations for the Lagrange multipliers A; has non-trivial solutions only if the determinant of the mx m 
matrix G is zero. The equation \Gjk \ = represents a condition on the parameter A. 



Now, 



jiven a solution to Eq. (A6), we take the expectation value {^C^)- Substituting from the right hand side of 



Eq. (A6), we find for this expectation value 



1^1 



A A 

1=1 1=1 



Xi\) 



1=1 

= A 



mo 



(AT) 
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In Eq. (A7) we have made use of the orthogonaUty of ^ to the Ws. We are also assuming that the function ^ is 
normaUzed. Thus, in solving for the value of A that satisfies Eq. (A6) we are also determining the effective values of 
the eigenvalues of the constrained problem. 



APPENDIX B: TRANSFER MATRIX 



When we are dealing with deep potential wells and large negative eigenvalues, the usual numerical integration of a 
differential equation over long intervals gives erroneous result. The reason is that as we integrate along a given path, 
the error starts growing exponentially and the longer is the distance, the more unreliable the final answer is. To avoid 
this, we only integrated numerically over half of a loop and with the help of the transfer matrix, we calculated the 
eigenf unctions in other regions. If we have a potential in the interval < x < L, and if the potential has reflection 
symmetry about x = L/2, we can express a solution ip{x) and its derivative at x = L/2 in terms of the two "primary" 
functions, $i(a;) and ^2{x), and their derivatives, at L/2, as follows 



^'(L/2) ; - $2(L/2) mL/2) ) [ ^'(0) ; 



Here, the functions ^i{x) and $2(2;) satisfy the equation in the interval. They also satisfy the boundary conditions 

$i(0) = 1 (B2) 

$'i(0) = (B3) 

$2(0) = (B4) 

$2(0) - 1 (B5) 

To find the function at the end of the interval, x — L, we reverse its sign in the middle and multiply by T~^, That is 
where 

R-(j_;) (B7) 

Effectively, we have a potential of the form of Figure || and an interval equal to 2L. 

To get anywhere along two loops, we use the appropriate combination of the above transfer and slope-reversing 
matrices. It is useful to construct a look-up table of the "fractional" transfer matrix t(.T), where 



where < x < L/2. With the use of this matrix, we can construct solutions throughout the interval. For instance, to 
obtain the solution in the interval L/2 < x < L, one makes use of the following relationship: 

^-Rt(L-.)RT-RTrSO\) (B9) 



Or, to obtain the solution in the interval 3L/2 < x < 2L, one utilizes: 



4^'ix) J - - x)RT-^RTRT-^RT ( J,^^^^ ) (BIO) 



Equation (BIO) may appear very complex. However, it saves a considerable amount of computational time and leads 
to a reliable answer. One good measure of accuracy of our answer is the Wronskian of the two independent solutions 
of the differential equation, which was found to be constant, as expected, along the interval. 
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FIG. 1. The Euler angles, 6, <t> and i). 




FIG. 3. The effective potential, V{s) in the operator in (3.1) 



de,,(s)/d& 




FIG . 4. The translational mode, d6ci{s) / ds. Also shown, as a dashed curve, is the effective potential, V{s), in the operator 
in (Ell). 



FIG. 5. The spectrum of the potential shown in Figure 




. L 0.1 . j 0.1 

FIG. 6. Alteration of the band structure as the result of a periodic array of constraints. The dotted curves are the bands in 
the absence of constraints. The solid curves represent the influence of constraints. 
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FIG. 7. Alteration of the band structure as the result of a periodic array of constraints, when there are two constraints per 
period. Here the solid curve is the result for \{q) when there are two constraints per period. The dashed curve represents the 
influence of one constraint, and the dotted curves are the bands in the absence of constraints. The solid curves represent the 
influence of constraints. 




FIG. 8. A schematic representation of a non-trivial supercoiled configuration of DNA stabilized by a protein armature. 




FIG. 9. Schematic of the nucleosome. 
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FIG. 10. The constrained spiral 
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FIG. 12. The portion of the compound spiral in Fig. Q that corresponds to a nucleosome in the 30 nm spiral. 




FIG. 13. The nuclesome configuration in this case. Note the the "histone" is in the form of a curved cylinder. The curvature 
is so as to fit into the two-loop structure that we investigate here. 



a) c) 




A= -6.04910501 X= -6.O491O207 

b) d) 




X= -6.04910207 x= -0.08066890 

FIG. 14. The four unstable eigenvalues of the fluctuation operator for the "nucleosome" configuration illustrated in Figure 
|l3|. The associated eigenvalues are shown immediately below the graphs. 
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FIG. 15. A figurative depiction of some of the constraints that were found to not stabilize the "nucleosome" configuration 
against fluctuations. 
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FIG. 16. A determinant associated with one of the constraints that does not stabilize the nuclesome configuration illustrated 
in Figure |l^. Note that this determinant passes through zero as a function of A for A < 0. Recall that zeros of the determinant 
arc proportional to energy eigenvalues of the fiuctuation operator. 




FIG. 17. A figurative version of the actual set of constraints that were utilized in this set of calculations. In this picture, the 
"histone" has been straightened out to resemble a the cylindrical shape that it actually takes. 
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\G(X)\ 




FIG. 18. The determinant yieldinK the eigenvalues of the cons trained fluctuation spectrum for the case of the constraints 
on the two loops in Figure la that are described in Section [VIII[ Note that this determinant as function of A does not pass 
through zero for any A < 0. Given that zeroes of the determinant are proportional to energy eigenvalues of the constrained 
fluctuation operator, we are assured that the nuclesome conflguration is stabilized against mechanical fluctuations. 
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